﻿ Cel mai bun estimator liniar nedeplasat, BLU De multe ori nu putem determina estimatorul MVU, chiar dacă el există Dacă, spre exemplu nu cunoaștem densitate de repartiție a datelor sau nu cunoaștem modelul de date, nu putem aplica metodele bazate pe CRLB sau pe statisticile suficiente Vom căuta un estimator ce depinde liniar de date, de forma 1N− ˆ = [] ∑ naxnθ 0n= Ponderile datelor se determină astfel încât estimatorul să fie nedeplasat și cu dispersia minimă Dacă estimatorul optimal e liniar ca funție de date, căutarea estimatorului BLU va conduce la estimatorul optimal Spre exemplu, pentru 2 ; 0,1, , 1, 0,xn A wn n N wnσ=+ = −∼N [][][] () estimatorul MVU este liniar funcție de date 1N− 1 ˆ [] ∑ MVUAx xn== 0nN= Dacă căutăm un estimator BLU pentru A, de forma 1N− ˆ [] ∑ nBLUAaxn= 0n= 1 vom găsi ponderile 1/N, adică chiar estimatorul MVU Nu întotdeauna estimatorul MVU este liniar după date De exemplu dacă se estimează componenta continuă a unui zgomot cu repartiție uniformă ; 0,1, , 1, 0, 2xn wn n N wnθ==−∼U () [][][] se arată că estimatorul MVU este 1N+ ˆ maxxnθ [] {} MVU= 2N Se poate lucra, în acest caz, cu estimatorul BLU suboptimal 1N− 1 ˆ [] ∑ BLUAxxn== 0nN= În unele probleme estimatorul BLU nu poate fi aplicat în mod direct Spre exemplu, pentru zgomotul alb, gaussian, de medie nulă, estimatorul MVU pentru dispersie este neliniar în date: 1^N− 221 xnσ= [] ∑ − 01nN= Dacă aplicăm datelor o transformare neliniară adecvată, care în cazul acesta este 2 ynxn= [][] se poate găsi estimatorul BLU 1^N− 2 ; 1 1ayn a Nσ==− 2[] nn∑ 0n= 1 Un alt exemplu în care se cere transformarea prealabilă a datelor este cel al repartiției log-normale, de parametru necunoscut θ 211⎧ ⎧⎫ exp ln , 0xn xnθ −−> [][] () ⎨⎬ ⎪ 22;xnpxnπ ⎩⎭ θ= [] [] () ⎨ 0 , xn⎪ 0 ≤ [] ⎩ Media repartiției log-normale se calculează cu 22 ∞∞ ⎧⎫⎧⎫ −− ln ln11xxθθ ()() ⎪⎪⎪⎪ =− exp expEx x dx dx=− {} ⎨⎬⎨⎬ ∫∫ 002222xππ ⎪⎪⎪⎪ ⎩⎭ ⎩⎭ 2 22yθ − () ∞∞ −−+ 22yyyθθ −− 2211y− Ex e e dy e dy== {} ∫∫ 22ππ −∞ −∞ 2⎡⎤ −+ 1yθ ⎢⎥ () ∞ 11⎢⎥ ⎣⎦ − θθ++ 2221y unde y=ln(x) sau x=eedye e== ∫ 2π −∞ Este evident că estimatorul nu poate satisface condiția de nedeplasare ˆ θ= {} BLUEθ În figură se prezintă densitatea de 3 repartiție log-normală pentru θ=0 2 Dacă facem transformarea prealabilă ln , 0xn xn> ⎧ [][] yn= [] ⎨ 0, 0xn≤ [] ⎩ densitatea de repartiție a variabilei aleatoare y se determină cu relația 2 ⎧⎫ − () 1yθ ⎪⎪ 2exp− ⎨⎬ yyθ − () 22eπ pxy− () () ⎩⎭ 21⎪⎪ == = () ypy e− 2edxπ dy ceeace înseamnă că y are o repartiție normală, de medie θși dispersie unitară ,1yθ∼N () Se poate aplica estimatorul pentru determinarea mediei unei variabile aleatoare normale 1N− 1 ˆ = 4[] ∑ BLUynθ 0nN= 2 Determinarea estimatorului BLU Condiţia ca estimatorul BLU să fie nedeplasat este 1N− ˆ EaExnθ, ==∀ [] {} {} ∑ nθθ 0n= Dispersia estimatorului se calculează cu 2 11NN−− ⎧⎫ ⎛⎞ ⎪ ˆ⎪ =− [] [] {} {}⎨⎬ ∑∑ nnDisp E axn aE xnθ ⎜⎟ 00nn== ⎝⎠ ⎪⎪ ⎩⎭ 2 1N− ⎧⎫ ⎡⎤ ⎪⎪ EaxnExn=− [] [] {} () ∑ n⎨⎬ ⎢⎥ 0n= ⎣⎦ ⎪⎪ ⎩⎭ Notăm cu a vectorul ce are ca și componente coeficienții estimatorului BLU T a aa a= ⎡⎤ 01 1N− ⎣⎦ 5 Estimatorul pentru parametrul necunoscut θ poate fi pus sub o formă compactă 0x⎡⎤ [] ⎢⎥ 1N− 1x [] ˆT⎢⎥ ax == = ⎡⎤ [] ∑ nNaa a axn− 01 1θ ⎣⎦ ⎢⎥ 0n= ⎢⎥ 1xN− [] ⎣⎦ Și dispersia se poate pune sub o formă compactă T ˆTT TT Disp E E Eθ=− −ax a x ax a x {}{} ()() {} {} T T EE E=− −ax x x x a {}{} ()() {} T T EE E=− −axxxxa {}{} ()() {} Matricea de covarianţă a datelor se determină cu relația de definiție T EE E=− −Cxxxx {}{} ()() {} 6 3 Dacă înlocuim matricea de covarianță dispersia se poate exprima cu ˆT Dispθ=aCa {} „Relaţia de constrângere” poate fi satisfacută doar dacă: Exn snθθ=∀, [][] {} Dacă condiția de mai sus nu este satisfăcută și, spre exemplu, avem θ Exn ne= [] {} condiţia de constrângere ne conduce la 1N− θ na eθθ=∀, n∑ 0n= relaţie imposibil de satisfăcut pentru orice valoare a parametruluiθ 7 Dacă datele au media dependentă liniar de parametrul necunoscut, adică Exn snθθ=∀, [][] {} atunci datele pot fi puse sub forma ⎛⎞ xnExn xnExn⎜⎟ =+− [] [][] [] {}{} ⎜⎟ de la medieabatere⎝⎠ sn wnθ =+ [] [] Privim astfel abaterea datelor de la medie ca pe un zgomot wn xn E xn=− [][][] {} Notăm cu s vectorul de date utile T 01 1ss sN⎡⎤=− s [] [] [ ] ⎣⎦ Condiția de constrângere devine 1N− asnθθθ=∀, [] n∑ 0n= sau, după simplificarea cu θ T as 81= 4 Pentru a determina estimatorul BLU, este necesar să minimizăm, după vectorul a, dispersia ˆT Dispθ=aCa {} cu satisfacerea condiţiei de constrângere T 1=as *** Trebuie să minimizăm o expresie, ținând seama de relația de constrângere Se construiește lagrangeanul TT 1Jλ=+−aCa as () Pentru o matrice simetrică, așa cum e și matricea de covarianță, avem regula de derivare T∂ 2=aCa Ca () ∂a Mai știm că T∂ TT∂ ==as sa s () ∂∂aa Putem acum deriva J J∂ 20λ=+=Ca s ∂a și găsim soluția 91− 2λ=−aCs () Se introduce soluția în relația de constrângere, din care se găsește λ 1TT− 12λ=−=as sC s () Avem deci 1λ −= 1T− 2sC s Substituim λ în expresia vectorului a șiobținem expresia vectorului ce minimizează dispersia, satisfăcând, în acelașitimp, șirelația de constrângere 1− Cs =a 01T− sC s Substituind vectorul determinat găsim minimul dispersiei estimatorului 11 1TT−− − 1sC CC s sC s T ˆ == = =aCa {} 00122TDispθ 11TT−−− sC s sC s sC s ()() Estimatorul BLU determinat este 1T− sC x ˆT ax θ== 10TBLU− 10sC s 5 Vom exemplifica cu un caz în care componenta continuă, A, necunoscută, este afectată de un zgomot alb, de medie nulă și dispersie (putere) cunoscută, dar căruia nu-i știm densitatea de probabilitate Deoarece ; [ ] [ ]; n=0,1, ,N-1Exn A xn A wn==+ [] {} vectorul de semnal se ia T 11 1==s1 [] Zgomotul w[n] fiind alb, matricea de covarianță are o formă diagonală 2 =CI uσ Aplicând formula stabilită pentru determinarea estimatorului BLU obținem 1 T 1Ix 21TNu− 11x ˆσ Axn x=== = [] ∑ 1T 11 0TnN= 2u1I1 σ că, oricare ar fi repartiția datelor, estimatorul BLU al 11ceeace înseamnă componentei continue este media eșantion În stabilirea relației anterioare am ținut seama de expresiile a două produse 0x⎡⎤ [] ⎢⎥ 1N− 1x [] T == 11 1xn⎢⎥1x [] [] ∑ ⎢⎥ 0n= ⎢⎥ 1xN− [] ⎣⎦ 1⎡⎤ ⎢⎥ 1 T 11 1N⎢⎥ == 11 [] ⎢⎥ ⎢⎥ 1 ⎣⎦ Dispersia estimatorului BLU se obține înlocuind în formulă Ea este 2 1σ ˆ Disp A== {} 1N T 2u1I1 σ A ar fi densitatea de probabilitate a zgomotului w[n] ce afectează12oricare 6 Vom considera acum acelaşi semnal continuu A, dar afectat de un zgomot de medie nulă, cu eșantioanele necorelate,dar în care dispersiile eşantioanelor de zgomot diferă între ele 2 Disp w n Nnσ=−= 0,1, , 1, [] {} n Deoarece eșantioanele de zgomot sunt necorelate, toate elementele matricei de covarianță sunt nule, cu excepția celor de pe diagonala principală care sunt dispersiile eșantioanelor Se dau matricea de covarianțăși inversa ei 22− ⎡⎤⎡ ⎤ 0000 0 0σσ ⎢⎥⎢ ⎥ 22− 111000 0σσ− ⎢⎥⎢ ⎥ CC == ⎢⎥⎢ ⎥ ⎢⎥⎢ ⎥ 22− −− 1100 0 0NNσσ ⎣⎦⎣ ⎦ 13 Pentru determinarea estimatorului BLU calculăm mai întâi numărătorul expresiei sale 2− 000xσ ⎤ ⎡⎤ [] 0⎡ ⎢⎥⎢⎥ 2− [] 11100Txσ− ⎥ ⎢⎥ 11 1⎢=1Cx [] ⎢⎥⎢⎥ ⎢⎥ ⎢⎥ 2− 100xNσ − [] 1N− ⎣⎦⎣⎦ 0x⎡⎤ [] 2⎢⎥ 0σ ⎢⎥ 1x⎢⎥ [] 1N− [] 2xn⎢⎥ == [] ∑ σ1211 1 ⎢⎥ 0nnσ= ⎢⎥ ⎢⎥ 1xN− [] ⎢⎥ 2⎢⎥ σ ⎣⎦ 141N− 7 și apoi numitorul 2− ⎡⎤ ⎡⎤ 0100σ ⎢⎥ ⎢⎥ 21N− − 111001Tσ− ⎢⎥ = 11 1⎢⎥ =1C 1 [] 2∑ ⎢⎥ ⎢⎥ 0nnσ= ⎢⎥ ⎢⎥ 2− ⎣⎦ 1100Nσ− ⎣⎦ Putem da acum expresia estimatorului BLU pentru A 1N− xn [] 2∑ 0nnσ ˆ= A= 1BLUN− 1 2∑ = 0nnσ 15 Se poate calcula, cu cele stabilite deja, dispersia estimatorului 11 ˆ == {} 21BLUNDisp A− 1S 2∑ 0nnσ= Coeficienții de ponderare a datelor, elemente ale vectorului a, sunt 11 2na= nSσ Estimatorul BLU pentru cazul parametrului vector necunoscut Dacă avem de estimat BLU p parametri, aceștia se organizează sub forma unui vector T ⎡⎤ θ = θ 112ppθθ × ⎣⎦ Expresiile celor p estimatori BLU sunt de forma 1N− ˆ == [] ∑ ,, 1, 2, ,iinaxn i pθ 16 0n= 8 Explicit, cele p relații de determinare a estimatorilor BLU, pentru cele p componente sunt ˆ =+++ − [][][] 11,0 1,1 1,101 1Nax ax a xNθ− ˆ =+++ − [] [] [ ] 22,0 2,1 2,101 1Nax ax a xNθ− ˆ =+++ − [] [] [ ] ,0 ,1 , 101 1ppp pNax ax a xNθ− Aceste p relații se pot rescrie și sub o formă matriceală ˆ ⎡⎤ ⎡⎤ ⎡⎤ [] 11,0 1,1 1, 10Naa axθ − ⎢⎥ ⎥ ⎢⎥ ˆ⎢ aa axθ [] ⎢⎥ 2,0 2,1 2, 121N− ⎢⎥ ⎢⎥ = ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ − ⎢⎥ [] − ,0 ,1 , 1ˆ1pp pNaa axNθ ⎣⎦ p⎣⎦ ⎣⎦ xA × 171ˆNpN× θ 1p× sau, cu notațiile matriceale obișnuite ˆ =θAx Pentru asigurarea condiției de nedeplasare a estimatorului este necesar să avem ˆ =∀θ Axθθ EE=, {} {} relație care poate fi îndeplinită pentru orice vector θ dacă vectorul mediei datelor are forma =xHθ {} 11NppNE×× × Substituind, se obține condiția ca să avem un estimator vector nedeplasat AHθ=θ conditie care afirmă că produsul AH este o matrice unitate =AH I pNNp upp××× Dacă vectorul linie din matricea A este T ⎡⎤ ==a ,0 ,1 , 1 1, 2, ,iii iNaa a i p− 18⎣⎦ 9 matricea A se poate pune sub forma T ⎡⎤ a ⎡⎤ 1,0 1,1 1, 11Naa a− ⎥ ⎢⎥ T⎢ a 2,0 2,1 2, 12Naa a− ⎢⎥ ⎢⎥ == A ⎢⎥ ⎢⎥ ⎥ ⎢⎥ T⎢ a ,0 ,1 , 1pp pNpaa a− ⎢⎥ ⎣⎦ ⎣⎦ În structura matricei H evidențiem vectorii coloană care o compun ⎡⎤ Hhh h 12p= ⎣⎦ Cu aceste notații, condiția de constrângere ia forma T ⎡⎤ a ⎤ 110 0⎡ T⎢⎥ ⎢⎥ 201 0a ⎢⎥ ⎢⎥ ⎡⎤ AH h h h 12p== ⎣⎦ ⎢⎥ ⎢⎥ ⎢⎥ T⎢⎥ 1900 1a ⎢⎥ p ⎣⎦ ⎣⎦ Putem deduce că avem 1, ij= T⎧ ==ah = ⎨ ,; , 1, 2, ,jijiij pδ 0, ij≠ ⎩ Conform cu formulele stabilite pentru parametri necunoscuți scalari, aplicate fiecărei componente, rezultă formele celor p estimatori BLU ˆT ; 1, 2, ,ipθ==ax ii și cele p dispersii ale estimatorilor 1T− ˆ ; 1, 2, ,Disp i pθ =aC a {} iii= Matricea parametrilor optimi se poate pune sub forma 1− − 11TT− AHCHHC () 0= vectorul estimatorilor BLU fiind 1− 11TT−− ˆ θHC H HC x () BLU= Matricea de covarianță a estimatorilor este 1− 1T− CHCH ) ˆ= 20( θ 10 Estimatorul BLU devine estimator MVU pentru modelul liniar generalizat, în care am admis că : =+xHθ ww 0,C∼N ; () media datelor fiid liniară în vectorul parametru θ, cerință pentru a putea stabili existența unui estimator BLU Din modelul de mai sus rezultă E=xHθ {} Teorema Gauss-Markov Dacă datele x urmează un model linia generalizat =+xHθ ww 0,C∼N ; () în care HNxpeste o matrice de observare cunoscută, θpx1 este un vector necunoscut, format cu paramerii ce trebuie determinațiiar wNx1 este un vector de zgomot de medie nulă și matrice de covarianță C, dar având o densitate de probabilitate oarecare, necunoscută, estimatorul BLU pentru vectorul θ este 1− 11TT−− ˆ θHC H HC x () BLU= și are matricea de covarianță C, dată de expresia 1− 1T− CHCH () ˆ= θ 21 Dispersiile componentelor vectorului de estimare sunt minime și au valorile 1− 1T− ⎤ ˆ⎡ HC Hθ () {} iDisp= ⎥ ⎣⎦ ,ii⎢ și sunt extrase de pe diagonala principală a matricei de covarianță a estimatorului vector *** SISTEMUL TDOA Vom analiza un exemplu privind localizarea unei surse ce emite semnale radio, care sunt recepţionate de mai multe antene fixe Semnalul radio este emis, de exemplu, de către un avion și este receptionat de trei antene fixe Calculând diferențele timpilor de sosire a undelor la antene, putem localiza sursa de semnal la intersecția a două hiperbole În figură se prezintă un caz particular, pentru care undele sosesc în același timp la cele trei antene Sursa se află la intersecția mediatoarelor segmentelor ce 22 unesc antenele 11 Vom considera un caz mai general, în care recepția se face cu N antene Dacă sursa emite un semnal la momentul (necunoscut) T0, recepț ia are loc, la antena i la momentul iR =++ = − 0 ; 0,1, , 1iitT i Nε c Am notat cu c viteza undelor electromagnetice, cu Ridistanța dela sursă i este cauzat de procedeul de măsurare, eșantioanele sale antena i; zgomotul ε având media nulă și dispersia cunoscută iar zgomotele sunt necorelete Vectorul parametrilor necunoscuți este constituit din coordonatele sursei de semnal T θ [] ssxy= Distanța de la sursă la antena i este 22 Rxx yy=−+−; ()( ) ssii i − 0,1, , 1iN= 23 În sistemele de urmărire (tracking) ale sursei radio dispunem, din măsurări efectuate anterior, de o poziţie “nominală” a sursei, apropiată de poziția ei efectivă Pentru “liniarizarea” problemei vom considera ca parametru vector T Δθ [] ssxy=Δ ale cărui componente au valori mici, reprezentând erorile geometrice de localizare 24 xxx y yyΔ= − Δ= −; ssn s sn 12 Pentru liniarizare vom dezvolta distanța sursă-antenă în serie Taylor, de ordinul unu − nniixx yy− 0,1, , 1;NRR x y i− Δ+ Δ = ssini≅+ ni niRR Dacă substituim această formă liniară a distanței în expresia timpilor obținem nnni i iRxx yy−− 0,1, , 1;NtT x y iε− ≅+ + Δ+ Δ+ = 0ssii ni niccR cR Din figură se vede că putem scrie relațiile nniixx yy−− cos0,1, , 1; sin ; Niαα− = ii== RR ii care introduse în expresia timpului ne conduc la cos sinRαα ni i i 0,1, , 1;NtT x y iε−≅+ + Δ+ Δ+ = 0ssii cc c 25 Notăm niR 0,1, , 1;Ntiτ−=− = ii c șiobținem cos sinαα ii 0,1, , 1;NTxyiτε−≅+ Δ+ Δ+ = 0ssii cc Cum momentul emiterii semnalului radio de către sursă nu se cunoaște, îl vom elimina, calculând N-1 diferențe, care pot fi măsurate ξττ ξ ττ ξ τ τ=− = − = −; ; ; 110 2 21 1 1 2NNN−−− Diferențele (măsurabile) sut dependente de erorile de determinare αα−− 11cos cos sin sinii iiαα−− =Δ+Δ+− ;xyξεε 1ssiii− cc = 1, , 1Ni− 26 13 Cele N-1 relații se explicitează ααα−− 10 10cos cos sin sinα ε=Δ+Δ+− 110ssxyξε cc αα−− 21 21cos cos sin sinαα ε=Δ+Δ+− 221ssxyξε cc αα−− −− −− 12 12cos cos sin sinNN NNαα ε=Δ+Δ+− −−− 112ssNNNxyξε cc și se trec sub formă matriceală ααα εε−− − ⎡⎤⎡ ⎤ ⎡ ⎤ 11010 10cos cos sin sinξα ⎢⎥⎢ ⎥ ⎢ ⎥ −− − ααα εε Δ ⎡⎤ 22121 21cos cos sin sinsxξα ⎢⎥⎢ ⎥ ⎢ ⎥ =+ ⎢⎥ ⎢ ⎥ ⎢ ⎥ Δ sy⎢⎥ ⎣⎦ ⎢⎥⎢ ⎥ ⎢ ⎥ θ αα εε −− − 21cos cos sin sinξα× −−−−− −− 11212 12NNNNN NNα ⎣⎦⎣ ⎦ ⎣ ⎦ wHξ 111211NNN−× −× −× () () ) ( Cu notații matriceale uzuale am stabilit relația 27 ξ=Hθ +w Se vede că ε− − ⎡⎤ ⎡⎤ ⎡⎤ 10 01100 00εε ⎢⎥ ⎢⎥ ⎢⎥ − ε − 21 10110 00εε ⎢⎥ ⎢⎥ ⎢⎥ == = wAε ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ − ε − ⎣⎦ −− − 12 10000 11NN Nεε ⎣⎦ ⎣⎦ wε A 111NN− 1NN−×−× ) ( ) ( Matricea de covarianță a zgomotului de măsurare este de forma 2 σ=CI uε iar a zgomotului rezultat în urma prelucrării datelor T TT T T EEE===CAεAεAεεAAεεA ()() {}{} {} 2TT σ ==AC A AA ε Aplicăm relațiile de determinare a estimatorului BLU pentru vectorul parametrilor necunoscuți T Δθ [] ssxy=Δ 28 14 Avem, pentru estimatorul BLU 1− 1111−− ˆTT TT⎡⎤ θ HAAHHAAξ ()() 22= ⎢⎥ σσ ⎣⎦ 1− 111−− 2TT T T ⎡⎤ = HAA H H AA ξ () () σ2 ⎢⎥ ⎣⎦ σ 1− 11−− TT TT = HAA H HAA ξ ⎡⎤ () () ⎢⎥ ⎣⎦ ce are matricea de covarianță 1− 1− 2TT ⎡⎤ σ=CHAAH () ˆ⎢ ⎥ θ ⎣⎦ Dispersiile fiecărui estimator scalar se calculeză cu elementele diagonalei principale a matricei de covarianță 1− ⎡⎤ 1− 2TT⎡⎤ ˆ =HAA H ⎢⎥ () {} iDispθσ ⎢⎥ ⎣⎦ ⎢⎥ ⎣⎦ 29,ii Vom analiza un caz particular, în care se utilizează N=3 antene, coliniare, ca în figură Luăm o poziție nominală particulară, indicată în figură Pentru valorile funcțiilor trigonometrice implicate în acest caz avem ααα== 00cos cos sin sinα == 11cos 0 sin 1αα =− = ααα 22cos cos sin sinα Rezultă forma particulară a matricei de observare 0cos 1sin cos 1sin11αααα−− −− ⎡⎤⎡⎤ ==H ⎢⎥⎢⎥ 30cos 0 sin 1 cos sin 1ccαα αα −− − − − ⎣⎦⎣⎦ 15 Matricea A are forma 110− ⎡⎤ A = ⎢⎥ 011− ⎣⎦ Calculăm produsul 10− ⎡⎤ 110 2 1−− ⎡⎤ T⎡⎤ ⎢⎥ =AA 11=− ⎢⎥ ⎢⎥ 011 12⎢⎥−− ⎣⎦ ⎣⎦ 01⎢⎥ ⎣⎦ determinantul produsului T 41 3=−=AA și matricea inversă a matricei produs 121− 1⎡ T⎤ AA = () ⎢⎥ 123 ⎣⎦ 31 Determinăm expresia produsului matriceal α−− − − 1cos cos 2 1 cos 1 sin1α− ⎡⎤⎡⎤⎡⎤ TTαα HAA H = () ⎢⎥⎢⎥ 2⎢⎥ 1 sin sin 1 1 2 cos sin 13cαα αα −−− − ⎣⎦⎣⎦⎣⎦ 3cos 3cos cos 1 sin1αα −−−− αα ⎡⎤⎡⎤ ⎢⎥⎢⎥ 2= 1 sin sin 1 cos sin 13cαα −−−− αα ⎣⎦⎣⎦ 2 ⎡⎤ 2cos 0α 1 ⎢⎥ 222= ⎢⎥ 01sincα − () 3⎢⎥ ⎣⎦ ce are determinantul 124− 2TT α=−HAA H () () 4cos 1 sinα 3c 32 16 Calculăm inversa produsului matriceal 22⎡ ⎤ − () 141sin 0α 2− ⎢⎥ 131c− 3TTc⎡⎤ HAA H = ⎥ () 22⎢ ⎢⎥ ⎣⎦ 224cos 1 sinαα ⎢⎥ − () 0cosα 2⎢ c⎥ ⎣⎦ 1⎡ ⎤ 0 2⎢ ⎥ 22cosα ⎥ = c⎢ 3 ⎢⎥ 0 2⎢ ⎥ 21 sinα − () ⎣⎦ 33 Calculăm acum factorul matriceal din estimatorul BLU, ce ponderează măsurările diferențelor 1− 11−− TT TT ⎡⎤ HAA H HAA () () ⎢⎥ ⎣⎦ 1⎡⎤ 20⎢⎥ αα −− ⎡⎤⎡ ⎤ 22coscos cos 2 111α = c⎢⎥ ⎥⎢ ⎥ 3⎢ −− 1sin sin 1 1 23cαα ⎢⎥ ⎣⎦⎣ ⎦ 20 ⎢⎥ − 21 sinα () ⎣⎦ 11⎡⎤ −− ⎢⎥ 212cos 2coscαα ⎡⎤ = ⎢⎥ 33123⎢⎥ ⎢⎥ ⎣⎦ − ⎢⎥ −− 21 sin 21 sinαα ()() ⎣⎦ 11⎡⎤ −− ⎢⎥ 2cos 2cosαα ⎢⎥ c= 11 ⎢⎥ − ⎢⎥ 21 sin 21 sinαα −− ()() ⎣⎦ 34 17 Se obține estimatorul BLU vector 11⎡⎤ −− ⎢⎥ Δ ⎡⎤ αα ⎡⎤ 12cos 2cossxξ c⎢⎥ = ⎢⎥ ⎢⎥ 11yξ ⎢⎥ ⎣⎦ ⎣⎦ 2sΔ − ⎢⎥ −− 21 sin 21 sinαα ()() ⎣⎦ din care se deduc estimatorii BLU scalari c ξΔ=− + xξ () 12s 2cosα c ξΔ= − () 12syξ 21 sinα − () Matricea de covarianță a estimatorului vector este 1⎡ ⎤ 210− ⎢⎥ 12cos− α 222TT⎡⎤ ⎢⎥ =CHAAHσσ () ˆc= 3⎢⎥ θ ⎢⎥ ⎣⎦ 20⎢ 3521 sin⎥ α − () ⎣⎦ De pe diagonala principală a matricei de covarianță se extrag dispersiile estimatorilor scalari 2 c Disp xΔ=; {} 2s 2cosα 22 3cσ Disp yΔ= {} 2s 21 sinα− () Dispersii mici se obțin (vezi figura) pentru 0α→ și ele sunt 2 c Disp xΔ=; {} 2s σ 22 3cσ {} sDisp yΔ= 2 36 18